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We describe two ideas useful in the construction of highly improved quark actions for simulations on coarse 
lattices: (i) Field transformations solve the doubler problem without destroying tree-level improvement for on- or 
off-shell quantities. The simplest example is the Sheikholeslami-Wohlert (clover) action. Going to the next order 
of improvement yields the class of D234 actions, (ii) Anisotropic lattices with a t < a s are useful because they push 
up the energy of unphysical branches of the dispersion relation (which are generic to highly improved actions), 
allow accurate mass determinations for particles with bad signal/noise properties (glueballs, P-state mesons), and 
enable one to simulate heavy quarks within a relativistic framework. We present first simulation results for the 
quenched light hadron and charmonium spectra obtained with a D234 action on anisotropic lattices. 



1. Introduction 

The last few years have seen a revival of the 
Symanzik improvement program [1-3]. The im- 
petus was largely provided by tadpole improve- 
ment (TI) [4] , which yields a simple first estimate 
of the large perturbative corrections to the tree- 
level coefficients in an improved action. In the 
case of non-relativistic QCD [5] and pure glue [6], 
the TI of tree-level coefficients has lead to accu- 
rate results on lattices as coarse as a — OA feu. 

Improved glue has classical 0(a A ) errors, and 
after TI any remaining quantum errors (of which 
the leading 0(a 2 a) piece is known [2] and can be 
eliminated) seem to be small [6]. For quarks, the 
Sheikholeslami-Wohlert (SW) action [3], which 
has classical 0(a 2 ) errors, is the most improved 
action that has been used in large-scale simula- 
tions. For full QCD simulations to be possible on 
coarse lattices, it therefore seems mandatory to 
design more highly improved quark actions [7,8]. 
Since exploratory simulations will have to be per- 
formed in the quenched approximation, for which 
the a = hadron spectrum is not a priori known, 
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it is important to have simple improvement tests 
for a quark action. As such we propose: 

(i) The "effective velocity of light" c(p) (defined 
by c(p) 2 = (E(p) 2 - £(0) 2 )/|p| 2 ) of various 
hadrons. Its deviation from 1 for non-zero masses 
and momenta is a measure of scaling violations. 

(ii) The "r-test". For actions with Wilson and 
clover terms one can obtain a rough estimate 
of the correctness of the coefficient of the clover 
term, by seeing whether the spectrum is invariant 
under (small) changes of the Wilson parameter r. 

(iii) The scaling of dimensionless quantities such 
as mpf/nip, m&/m p , and J [9]. 

(iv) The lattice spacing dependence of the rho 
mass in units of one of the standard scales (string 
tension, charmonium S-P splitting, etc). This is 
a sensitive probe of scaling violations, but it re- 
quires that the scaling, systematic and measure- 
ment errors of the chosen scale be considerably 
smaller then those of the rho mass. 

Last year we presented results for the "isotropic 
D234" quark action [7], which has classical 0(a 3 ) 
errors. It gave very impressive results for all 
hadron mass ratios and the dispersion relation of 
mesons, much better than results from the SW ac- 
tion. However, the absolute value of the rho mass 
was almost identical to that of the SW action [10] 
for the three lattice spacings where they could be 
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compared. To understand this somewhat incon- 
gruous behavior we have started to explore a va- 
riety of improved quark actions. Our aim is to 
disentangle the various sources of scaling errors. 

One possible source of scaling violations is the 
unphysical branches in the quark dispersion rela- 
tion that are generically present in actions im- 
proved beyond 0(a). This has lead us to the 
study of actions on anisotropic lattices, and, in 
particular, to devise a simple procedure to con- 
struct actions that are tree-level improved to any 
order and doubler-free on an arbitrary lattice [11]. 
Working on anisotropic lattices with a t < a s has 
several advantages: 

(a) By decreasing a t one can decouple the un- 
physical branches by pushing up their energy. 

(b) Larger and more easily identifiable effective 
mass plateaux lead to more accurate and confi- 
dent mass determinations. This is important for 
particles with bad signal/noise properties, like P- 
state mesons and especially [12] glueballs. 

(c) Small at allows one to simulate heavy quarks 
within a relativistic framework without the pro- 
hibitive cost of a fine spatial lattice. 

The disadvantage of anisotropic lattices is that 
there are more independent coefficients in an im- 
proved action, and they have to be tuned to re- 
store space-time exchange symmetry. One hopes 
that suitable TI gives reasonable estimates of 
(most of) the coefficients, also on anisotropic 
lattices, and our tests indicate that in general 
([7,13], Sect. 3) this appears to be the case, ft 
is always better, of course, to tune them non- 
perturbatively. We briefly discuss this possibility 
in sect. 4. Tuning the clover term is important 
already for the isotropic case. It could solve the 
one problem of the isotropic D234 action, its low 
rho mass [7]. 

2. Classically Improved Quark Actions 

There is a simple procedure to construct tree-level 
improved quark actions on an arbitrary lattice 
without a doubler problem. Let a M be the lattice 
spacings of a generic hypercubic lattice, [do = a t ; 
when the spatial lattice spacings a, are identical 
they will be denoted by a s , such a lattice with 
£ = a s /a t will be referred to as a '£:1 lattice']. 



The idea is described in four steps: 

(i) Start with a naive (improved) fcrmion action, 
$ c M c ip c , M c = f c + m c , where V CM is an (im- 
proved) discretization of the continuum Dirac op- 
erator Dp, differing from it at C(a"), say. The 
subscript 'c' stands for 'continuum- like', since this 
action is manifestly improved to 0(a n ). How- 
ever, this action will have doublers. The simplest 
choice of V CM , which leads to the SW action, is 
the standard anti-hermitean covariant derivative, 

V^(x)^^-[^ + ,-^]=D, + 0(al) (1) 

with = U±^(x)ip(x ± fi). We will be inter- 
ested in higher orders of improvement, as exem- 
plified by the improved derivative 

V c „ = V M - 1 al V, A, = D, + . (2) 

Here A p is the standard second order derivative, 

A^(x) ee ^ + 1>-» - 2 V>(*)] • (3) 
a fj, 

(ii) To cast the action in a form where doublers 
can easily be eliminated, perform a field transfor- 
mation ip c = Clip, tp c = tpCl, so that ipcMcipc = 
•ipMnip, where Mq ee ClM c Q. We choose f2 = 
(when acting to the right) and 

fl SI = 1 - ±ra (% - m c ) + 0(a n ) , (4) 

in terms of the (initially) free parameter r, the 
Wilson parameter. The transformed fermion op- 
erator Mq reads at this point 

M n = m c (l + ±ra m c ) + f c -Ka f 2 c +O(a n ) .(5) 

This action still has doublers. 

(iii) To remove the doublers, use 

y c 2 = A,-^alAl+ l -a-F+... (6) 

to the appropriate order to eliminate V M in terms 
of A^. Here a ■ F = ^^cr^F^ is the clover 
term, with F^ v a lattice representation of the 
field strength, improved to the appropriate order. 
Since in momentum space A M does not vanish at 
the edge of the Brillouin zone, the action obtained 
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by the above truncation is doubler-free. This ac- 
tion, denoted by tpMip, is on-shell improved up to 
0(a n ) errors. If one is only interested in spectral 
quantities, one can stop here, 
(iv) To (classically) also improve non-spectral 
quantities, undo the field transformation, which, 
e.g., amounts to using the propagator G = 
£lM~ 1 il [15]. Note that this step does not rein- 
troduce the doublers, since it was the truncation 
in step (iii), not the field transformation per se, 
that eliminated them. 

We emphasize that in this approach one does 
not have to separately check the improvement of 
the interactions; it follows from the fact that they 
are improved for the naive action we started with, 
and this fact is not affected by the field transfor- 
mation (the Jacobian of the latter can be ignored 
to the order we are working). 

Let us explicitly write down the actions M ob- 
tained for n = 2 and 4: The SW action 

M sw = m a + f -Kaoly^A^ + ^a-Fj , (7) 

and the D234 action 

M D234 = m + 7„V M (1 - & M ajA„) 

+0^1 -i r ao(E M A^ + ia-F), ( 

where m = m c (l + ^ra m c ) in both cases, and 

^ = 1/6, c p = rao/(24a M ) . (9) 

We now choose r to lead to as few and high 
lying unphysical branches in the free dispersion 
relation as possible. The general D234 action (8) 
has four branches E = E(p) (not counting the 
particle anti-particle symmetry E —E), of 
which three are unphysical. One easily shows: 

• For b = 2c there will be at most two unphysi- 
cal branches. 

• If furthermore r = 1 — 2&o or bo = there is (at 
most) one unphysical branch. 

• There are no unphysical branches if and only if 
r = l,& = c o = 0. 

For the SW action r = 1 is therefore the canon- 
ical choice. For a D234 action obtained as above 
one will have at least one unphysical branch; ex- 
actly one if one chooses r = | and c = in 
addition to bo = |. To obtain an isotropic D234 




Figure 1. a s E(p) vs a s \p\, with p cx (1,1,0), for 
masslcss SW action on a 1:1 lattice (solid), and 
D234(|) on a 2:1 lattice (dashed). Continuum 
fermions (thin solid) are shown for comparison. 



action (on an isotropic lattice) one must choose 
all b^ — i, c M = and r = |. Since the c M vio- 
late (9), this action [7] has 0(a 3 ) errors. Like the 
r = 1 SW action, it can be coded very efficiently 
using the "projection trick" . 

Requiring just one unphysical branch is not 
compatible with (9). Suitably modifying c 
(and/or b at O(a )) to have just one branch leads 
to actions with O(ao,aj l ) errors. On anisotropic 
lattices with a s /a t > 2, say, one can presum- 
ably tolerate 0{a\) errors. In the free case it is 
even possible to avoid these 0{a\) errors; this is 
achieved [11] by a more complicated field trans- 
formation than (4). It leads to a D234 action 
where (some of) the coefficients have a slight mass 
dependence. Its spatial coefficients satisfy (9) and 
for the temporal ones, r = 1 — 2&o = 1 — 4co = 
|(1 + |aom c )/(l + j^aom c ). We will denote this 
action as D234(|), labelling it by the value of r 
at m c = 0. This is the action used in our prelimi- 
nary simulations. For details and some variations 
of the above procedure we refer to [11]. 

In figures 1 and 2 we compare the free quark 
dispersion relation of D234(|) on a 2:1 with that 
of SW on a 1:1 lattice. Note particularly the (ex- 
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u -° t i i i i i lJ 

o l i 2 3 4 0.2 0.4 0.6 0.8 1.0 1.2 

|p|/GeV 

Figure 2. As in figure 1, for masses such that 

a s E(0) = 1 in all cases. Figure 3. c(p) at various momenta for the 

D234(|) action on a 2:1 (pion •, rho o) and the 
SW action on a 1:1 lattice (pion o, rho x). For 
both actions m p /m^ « 1.47 and a s « 0.32 fm. 

pected) dramatic improvement observed in the 
heavy mass case. 



3. Simulation Results 

The anisotropic Symanzik-improved gluon ac- 
tion as well as the TI prescription we used is de- 
scribed in [12]. In addition to 0(af) and 0(a 2 8 a) 
errors, our gluon action has O(af) errors from 
leaving out terms that give rise to unphysical 
branches. For anisotropies a s /a t > 2 the O(af) 
errors seem to be small. However, we have no- 
ticed that, as the coupling is varied, the ratio of 
string-tension to charmonium determinations of 
a s shows variations of about 5%. This may be due 
to small violations of rotational invariance, and 
we are in the process of non-perturbatively tun- 
ing the gluon action to further reduce them [13]. 

For our preliminary investigation of quenched 
QCD on anisotropic lattices, we have chosen to 
study the D234(|) action (with the same TI pre- 
scription as for the glue) on 2:1 lattices for light 
quarks, and on 3:1 lattices for heavy masses. We 
use the SW action in various comparisons; both 
our as well as the results of [10] use TI. For our 
general simulation methodology we refer to [7] . 

For light quarks we performed simulations at 
three couplings, /3 = 2.0,2.15,2.3, corresponding 



to lattice spacings a s = 0.41, 0.36, 0.32 fm, if one 
uses charmonium to set the scale [13]. The spatial 
size of our lattices is about 2.4 fm (smaller for 
large quark masses). 

The first question is whether there are large 
0(a ) effects violating space-time symmetry. As 
one can see in figure 3, after TI such effects arc 
very small. In addition, figures 3 and 5 show that 
D234(|) is also more continuum-like at momenta 
up to 1 GeV, indicating that with TI the C(a 2 ) 
terms have roughly the right coefficients. 

In figure 4 we see that the D234(|) rho mass, 
converted to physical units using the charmonium 
S— P splitting, is closer to the continuum value on 
coarse lattices than for SW. 

The value of J [9] in table 1 is around 0.39(1) 
already on the coarsest lattices (as for isotropic 
D234 [7]); SW reaches this value only for much 
finer lattices [10]). The m^lnip ratio appears 
rather high, especially at (3 = 2.3, where it is 
larger than that of [7] as well as [10]. However, 
this might be due to chiral fitting uncertainties 
and problems with the glue. 

Dramatic improvements are seen for the char- 
monium spectrum on anisotropic lattices. In fig- 
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Figure 4. The rho mass for the 2:1 D234(§) and 
1:1 SW [10] actions. 
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Table 1 

Comparison of 2:1 D234(f ) and 1:1 SW [10] 



Figure 5. Effective velocity of light c(p) for the 
D234(|) and SW actions on various lattices at 
the charm mass, |p| « 670 MeV. Each quadruplet 
represents, from left to right, the r\ c ("pion"), J/ip 
("rho"), "nucleon" and "delta". For clarity we 
have slightly displaced the x-coordinates around 
the actual a,. 
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Figure 6. Quenched charmonium hyper- fine split- 
ting (in MeV): D234(f ) 3:1 (•), SW 3:1 (o), SW 
1:1 (*), FNAL (o) [16], NRQCD (□) [5], experi- 
ment (x). 



ures 5 and 6 we show results for D234(|) on 3:1 
lattices and compare them with SW on 3:1 and 
1:1 lattices. As expected, isotropic SW shows 
0(1) scaling violations; SW on 3:1 is much better, 
and D234(§) on 3:1 is much better still. 

Overall, anisotropic lattices are clearly well 
suited to coarse-lattice studies of heavy quarks 
and glueballs [12]. We have also seen some im- 
provement in the rho mass for light quarks. How- 
ever, our present results have systematic uncer- 
tainties in the glue and scale setting, so that it 
remains to be seen how big such improvements 
will ultimately be. For improved quark actions 
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to be accurate on isotropic or anisotropic coarse 
lattices it might be necessary to eliminate 0(a) 
quantum errors by non-perturbative tuning of the 
coefficients. 

4. Non-perturbative Tuning and the r-Test 

Our simulations so far rely on Tl to estimate 
the coefficients in an action. We are now start- 
ing work on non-perturbative tuning. For an 
anisotropic action such a tuning might be nec- 
essary already at O(a ), since the spatial and 
temporal first order derivatives might renormalize 
differently. As figure 3 shows, after TI this effect 
is small; the remainder could easily be tuned by 
introducing a "bare velocity of light" into the ac- 
tion and demanding c(p) = 1 for small masses and 
momenta. 

At O(a) one has to tune the clover term. Let us 
for the moment restrict ourselves to the isotropic 
case, where there are only two terms at 0(a), 
which we write as -fKE M + f cr-F). The 
coefficient r, say, can be adjusted to any desired 
value by a field transformation. To eliminate 
quantum errors at 0{a) one therefore has only to 
tune the (relative) clover coefficient u). A tuning- 
method was presented in [14] for the SW action 
on Wilson glue, for which it was found that u (at 
zero quark mass) is significantly larger than the 
TI estimate. Note that if this were also the case 
for the D234 action, it might solve the problem 
of its low rho mass. 

An alternative method could be based on tun- 
ing u> so that the spectrum of the action is invari- 
ant under (small) changes of the Wilson r param- 
eter, cf. sect. 1. If u; were independent of r (as it 
is classically) this will only be true for the correct 
non-perturbative w. The fact that w presumably 
depends on r at the quantum level, complicates 
the conversion of the r-test from a rough con- 
sistency check into a full tuning method. How- 
ever, this is not a problem of principle, and we 
are presently exploring this idea. 

On anisotropic lattices there is the additional 
complication, for both tuning and the r-test, that 
the temporal and spatial parts of the O(a) terms 
can suffer a relative renormalization. Assuming 
that TI makes this renormalization small, as we 



have seen for the O(a ) terms, we can however 
also apply the r-test to anisotropic lattices. We 
have obtained preliminary results for the D234 
actions of sect. 3. On a 2:1 lattice with a s = 

0. 36 fm we measured the difference between the 
rho masses for r = | and r = 1 for (i) u> = 0, and 
(ii) the tree-level tadpole improved w. Tuning 
nip/m^ to be 1.3 in all cases, we find the change 
in the rho mass (in lattice units) to be 0.104(10) 
in the first case, and 0.012(4) in the second. This 
suggests that the correct ui is somewhat, but not 
dramatically larger than the TI estimate (even 
taking into account other results indicating that 
w increases with decreasing mass). 
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